Effective one body Hamiltonian of two spinning black-holes with 
next-to-next-to-leading order spin-orbit coupling 
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Building on the recently computed next-to-next-to-leading order (NNLO) post-Newtonian (PN) 
spin-orbit Hamiltonian for spinning binaries [1] we improve the effective-one-body (EOB) description 
of the dynamics of two spinning black-holes by including NNLO effects in the spin-orbit interaction. 
The calculation that is presented extends to NNLO the next-to-leading order (NLO) spin-orbit 
Hamiltonian computed in Ref. [2]. The present EOB Hamiltonian reproduces the spin-orbit coupling 
through NNLO in the test-particle limit case. In addition, in the case of spins parallel or antiparallel 
to the orbital angular momentum, when circular orbits exist, we find that the inclusion of NNLO 
spin-orbit terms moderates the effect of the NLO spin-orbit coupling. 

PACS numbers: 04.25.-g,04.25.Nx 



I. INTRODUCTION 

Coalescing black-hole binaries are among the most 
promising gravitational wave (GW) sources for the cur- 
rently operating network of ground-based interferomet- 
ric GW detectors. Since the spin-orbit interaction can 
increase the binding energy of the last stable orbit, and 
thereby leading to large GW emission, it is reasonable 
to think that the first detections will concern binary sys- 
tems made of spinning binaries. For this reason, there is 
a urgent need of template waveforms accurately describ- 
ing the GW emission from coalescing spinning black-hole 
binaries. These template waveforms will be functions of 
at least eight intrinsic real parameters: the two masses 
mi and 1112 and the two spin- vectors Si and S2. Be- 
cause of the multidimensionality of the parameter space, 
it seems unlikely for state-of-the-art numerical simula- 
tions to densely sample this parameter space. This gives 
a boost to develop analytical methods for computing the 
needed, densely spaced, bank of accurate template wave- 
forms. Among the existing analytical methods for com- 
puting the motion and the dynamics of black hole (and 
neutron star) binaries, the most complete and the most 
promising is the effective-one-body approach (EOB) [3- 
8]. Several recent works have shown the possibility of 
getting an excellent agreement between the EOB ana- 
lytical waveforms and the outcome of numerical simula- 
tions of coalescing black-hole (and inspiralling neutron- 
star [9, 10]) binaries. A considerable part of the current 
literature deals with nonspinning black- hole systems [11- 
18], with different (though not extreme) mass ratios (see 
in particular [19, 20]) or in the (circularized) extreme- 
mass-ratio limit [21-25] (notably including spin [26]). 

The work at the interface between numerical relativity 
and the analytical EOB description of spinning binaries 
has been developing fast in recent years. The first EOB 
Hamiltonian which included spin effects was conceived in 
Ref. [6] . It was shown there that one could map the 3PN 
dynamics, together with the leading-order (LO) spin- 
orbit and spin-spin dynamical effects of a binary systems, 
onto an effective test-particle moving in a Kerr-type met- 



ric, together with an additional spin-orbit interaction. 
In Ref. [27] the use of the nonspinning EOB Hamilto- 
nian augmented with PN-type spin-orbit and spin-spin 
terms allowed to carry out the first (and up to now, only) 
analytical exploratory study of the dynamics and wave- 
forms from coalescing spinning binaries with precessing 
spins. Recently, Ref. [2], building upon the PN-expanded 
Hamiltonian of [28], extended the EOB approach of [6] 
so to include the next-to-leading-order (NLO) spin-orbit 
couplings (see also Refs. [29, 30] for a derivation of these 
couplings in the harmonic-coordinates equations of mo- 
tion and Ref. [31] for a derivation using an effective field 
theory approach). Using this model (with the addition 
of EOB-resummed radiation reaction force [7, 22, 32]), 
Ref. [33] performed the first comparison with numerical- 
relativity simulations of nonprecessing, spinning, equal- 
mass, black-holes binaries. Then, building on Ref. [2, 6] 
and Ref. [34], Ref. [35] worked out an improved Hamil- 
tonian for spinning black-hole binaries. 

Recently, Hartung and Steinhoff [1] have computed the 
PN-expanded spin-orbit Hamiltonian at next-to-next-to- 
lcading order (NNLO), pushing one PN order further 
the previous computation of Damour, Jaranowski and 
Schafcr [28]. The result of Ref. [1] completes the knowl- 
edge of the PN Hamiltonian for binary spinning black- 
holes up to and including 3.5PN. 

This paper belongs to the lineage of Refs. [2, 6] and 
it aims at exploiting the PN-expanded Hamiltonian of 
Ref. [1] so as to obtain the NNLO-accurate spin-orbit 
interaction as it enters the EOB formalism. Note that, 
by contrast to Refs. [35] and [2], we shall not discuss 
here spin-spin interactions, nor shall we try to propose a 
specific way to incorporate our NNLO spin-orbit results 
into some complete, resummed EOB Hamiltonian. Al- 
though the Hamiltonian that we shall discuss here does 
not resum all the spin-orbit terms entering the formal 
"spinning test-particle limit" , we shall check that it con- 
sistently reproduces the "spinning test-particle" results 
of Ref. [35]. 

The paper is organized as follows: in Sec. II we re- 
call the structure of the PN-expanded spin-orbit Hamil- 
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tonian (in Arnowitt-Dcscr-Misncr (ADM) coordinates) of 
Ref. [1] and then we express it in the center of mass frame. 
Section III explicitly performs the canonical transforma- 
tion from ADM coordinates to EOB coordinates and fi- 
nally computes the effective Hamiltonian, and, in partic- 
ular, the effective gyro-gravitomagnetic ratios. In Sec. IV 
we discuss the case of circular equatorial orbits, we derive 
the test-mass limit and we exploit the gauge freedom to 
simplify the expression of the final Hamiltonian. 

We adopt the notation of [2] and we use the letters 
a, b = 1,2 as particle labels. Then, m a , x a = (x l a ), p a = 
(j>ai), and S a = (S a i) denote, respectively, the mass, the 
position vector, the linear momentum vector, and the 
spin vector of the ath body; for a / d we also define 
Yab = x a - x&, Tab = \r a b\, n a b = r ab /r ab , I • I stands here 
for the Euclidean length of a 3- vector. 

II. PN-EXPANDED HAMILTONIAN IN ADM 
COORDINATES 

We closely follow the procedure of Ref. [2] . The start- 
ing point of the calculation is the PN-expandcnd two- 
body Hamiltonian H which can be decomposed as the 
sum of an orbital part, H D , a spin-orbit part, H so (linear 
in the spins) and a spin-spin term H ss (quadratic in the 
spins), that we quote here for completeness but that we 
are not going to discuss in the paper. It reads 

#(x Q ,p a ,S a ) = if (x Q ,Pa) +ffso(x a ,Pa,S Q ) 



+ tf ss (x a ,p a ,S Q ). (1) 

The orbital Hamiltonian H a includes the rest-mass con- 
tribution and is explicitly known (in ADM-like coordi- 
nates) up to the 3PN order [36, 37]. It has the structure 

iJ (x a ,p a ) = ^m a c 2 +iJ oN (x Q ,p a ) 

a 

+ ~n -^olPN (X a , p a ) + —7 H o2 PN (x a , Pa ) 
C l C 4 

+ ^ifo3PN(x 0)Pa ) + O^V (2) 
The spin-orbit Hamiltonian H so can be written as 

H so (x a , p a , S a ) = ^ r2 a (xb, Pb ) • S a . (3) 

a 

Here, the quantity £l a is the sum of three contributions: 
the LO (oc 1/c 2 ), the NLO (oc 1/c 4 ), and the NNLO one 
(oc 1/c 6 ), 

o a (x h , P6 ) = ni (xi, P6 )+^ LO (x hj Pi ,)+n* NLO (x&, Pb ). 

(4) 

The 3-vectors £1^° and fi^ LO were explicitly computed 
in Ref. [28], while fl^ NLO can be read off Eq.(5) of 
Ref. [1]. We write them here explicitly for completeness. 
For the particle label a = 1, we have 
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and n^ NLO can be ob- 



The expressions for ^l 2 ° , 
taincd from the above formulas by exchanging the parti- 
cle labels 1 and 2. 

Let us consider now the dynamics of the relative mo- 
tion of the two body system in the center of mass frame, 
which is defined by setting P i + p 2 = 0. Following [2], 
we rescale the phase-space variables R = xi x 2 and 
P = pi = p 2 of the relative motion as 

R P P i p 2 

r = T^lvJ' p= — = — = , (6 

GM /i /i fx 

where M = mi + m 2 and fi = mim 2 /M. In addition, we 
rescale the original time variable T and any part of the 
Hamiltonian as 
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NR 



(7) 



where i? NR = H — Mc 2 denotes the "nonrelativistic" 
Hamiltonian, i.e. the Hamiltonian withouth the rest- 
mass contribution. As in [2] we work with the following 
two, basic combinations of the spin vectors: 



S = Si + S 2 = micai + ?n 2 ca 2 , 



m 2 
mi 



Si 



mi , 
m 2 



m 2 cai + mica 2 , 



(8) 
(9) 



where we have also introduced the Kerr parameters of 
the individual black-holes, ai = Si/(mic) and a 2 = 
S 2 /(m 2 c). We recall that in the formal 1 "spinning test 
mass limit" where, for example, m 2 — > and S 2 — > 0, 
while keeping a 2 = S 2 /(m 2 c) fixed, one has a "back- 
ground mass" M ~ mi, a "background spin" Sbckgd = 
A/cabckgd — Si = micai, a "test mass" [i ~ m 2 , and 
a "test spin" St cs t = S 2 = m 2 ca 2 ~ /ica tes t [with 
atest = Stcst /(a* c )]- Then, in this limit the combina- 
tion S ~ Si = ?7iicai ~ Afcabckgd = Sbckgd mea- 
sures the background spin, while the other combination, 
S* ~ mica 2 ~ Mca t cst = -^fStest//-* measures the (spe- 
cific) test spin a tes t = Stest /(/•*<?) . Finally, since the use 
of the rescaled variables corresponds to a rescaling of the 
action by a factor l/(GMfi), it is also natural to work 
with the corresponding rescaled variables 



ax _ 



for any label X (X= 1,2, ,*). 



(10) 



As noted in Ref. [2] this formal limit is not relevant for the physi- 
cally most important case of binary black holes, for which a.2 — > 
and ni2 — > 0. 
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Using the definitions (6)-(10), the center-of-mass spin- 
orbit Hamiltonian (divided by /x) in terms of the rescaled 
variables has the structure 

g so (r,p,S,S*)= gso(r ' P ' S ' S * } (11) 

M 

= iff£°o(r,P,S,S*) 

cr 

+ ^NLo(r,P,S,S*) 

+ ^NNLo(r,P,S,S*) + 0^, 

(12) 

and it can be written as 

iT so (r,P,S,S*) - ^ (gt DM (S,n,p)+g§P M (S*,n,p)) , 

(13) 



with the following definitions: v = fi/M is the sym- 
metric mass ratios and ranges from (test-body limit) 
to 1/4 (equal-mass case); the notation (Vi,'V^,'V£) = 
Vi • (V 2 x V 3 ) = eijkVtVjVf stands for the Euclidean 
mixed products of 3- vectors; n = r/|r|; gg DM and g§P M 
are the two (dimensionless) gyro-gravitomagnetic ratios 
as introduced (up to NLO accuracy) in [2]. These two 
coefficients parametrize the coupling between the spin 
vectors and the apparent gravito-magnetic field seen in 
the rest-frame of a moving particle. Their explicit ex- 
pressions including the NNLO contribution read 
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(14b) 



The label "ADM" on the gyro-gravitomagnetic ra- 
tios (14) is a reminder that, although the LO values are 
coordinate independent, both the NLO and NNLO con- 
tributions to these ratios actually depend on the defi- 
nition of the phase-space variables (r, p). In the next 
Section we shall introduce the two, related, effective 
gyro-gravitomagnetic ratios that enter the effective EOB 
Hamiltonian, written in effective (or EOB) coordinates, 
according to the prescriptions of [2] . 

III. EFFECTIVE HAMILTONIAN AND 
EFFECTIVE GYRO-GRAVITOMAGNETIC 
RATIOS 

Following Ref. [2], two operations have to be per- 
formed on the Hamiltonian written in the center of mass 



frame so to cast it in a form that can be resummed in 
a way compatible to previous EOB work. First of all, 
one needs to transform the (ADM) phase-space coordi- 
nates (x a , p a , S a ) by a canonical transformation compat- 
ible with the one used in previous EOB work. Second, 
one needs to compute the effective Hamiltonian corre- 
sponding to the canonically transformed real Hamilto- 
nian. Following the same procedure adopted in [2], we 
start by performing the purely orbital canonical trans- 
formation which was found to be needed to go from the 
ADM coordinates used in the PN-expandcd Hamiltonian 
to the coordinates used in the EOB dynamics. Since in 
Ref. [2] one was concerned only with the NLO spin-orbit 
interaction, it was enough to consider the lPN-accurate 
transformation. In the present study, because one is 
working at NNLO in the spin-orbit interaction, one needs 
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to take into account the complete 2PN-accurate canon- 
ical transformation introduced in [3]. The transforma- 
tion changes the ADM phase-space variables (r, p, S, S*) 
to (r', p', S, S*) and it is explicitly given by Eqs. (6.22)- 
(6.23) of [3]. To our purpose, we actually need to use 



the inverse relations r = r(r',p') and p = p(r',p'), so 
to replace (r, p) with (r',p') in Eq. (13). The needed 
transformation is easily found by solving, by iteration, 
Eqs. (6.22)-(6.23) of [3], and we explicitly quote it here 
for future convenience. It reads 
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As pointed out in [3], in the test-mass limit (y — > 0) 
one has r' J = [l + l/(2c 2 r)] r l , which is the relation be- 
tween Schwarzschild (r') and isotropic (r) coordinates 
in a Schwarzschild spacetime 2 . When this transforma- 
tion is applied to to the spin-orbit Hamiltonian in ADM 



coordinates, Eq. (13), one gets a transformed Hamilto- 
nian of the form H'(r', p', S, S*) = H' {v' , p', S, S*) + 
H' so (r',p',S,S*), with the NNLO spin-orbit contribu- 
tion that explicitly reads 



i^ LO (r',p',S,S*) = — US*,n',p' 
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(17) 



where we introduced the radial unit vector n' = r'/|r'|. With this result in hands, we can further perform on 

it a secondary purely spin- dependent, canonical transfer- 
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mation that affects both the NLO and NNLO spin orbit 
terms. This transformation can be thought as a gauge 
transformation related to the arbitrariness in choosing 
a spin-supplementary condition and in defining a local 
frame to measure the spin vectors. Such gauge condition 
can then be conveniently chosen so to simplify the spin- 
orbit Hamiltonian. This procedure was pushed forward, 
at NLO accuracy in Ref. [2]. In that case, the canonical 
transformation was defined by means of a 2PN-accurate 
generating function, that was chosen proportional to the 
spins and with two arbitrary (z^- dependent) dimension- 
less coefficients a(v) and b{v). Using rescaled variables, 
the NLO generating function of [2] reads 

G s2 pn = ^^P- {a{v){S,n',p') + b{v){S*,n', P ')) . 

(18) 

In Ref. [2] the parameters a(v) and b(v) were selected 
so to remove the terms proportional to p 2 in the final 
(effective) Hamiltonian. Let us recall that, at linear order 
in the G S 2pn, that was enough for the NLO case, the new 
Hamiltonian was computed as H" s °(y") = H' s °{y") — 
{H', G S 2pn}(u"), were we address collectively with y" = 
(r", p", S", S"*) the new phase space- variables. 

We wish now to introduce a more general gauge trans- 
formation such to act also on the NNLO terms of the 
Hamiltonian. To do so, in addition to the NLO part 
G S 2pn of the spin-dependent generating function men- 
tioned above, one also needs to introduce a NNLO con- 
tribution of the form 



G S 3PN — 
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(19) 



with six, arbitrary, independent dimensionless coeffi- 
cients. We shall then consider the effect of a spin- 
dependent generating function of the form G s = G S 2PN + 
G S 3PN- Since G 8 starts at 2PN order, it turns out that 
possible quadratic terms in the generating function are 
of order c -8 , i.e. at 4PN and thus are of higher order 
than the NNLO accuracy that we are currently consid- 
ering in the spin-orbit Hamiltonian. The consequence is 
that the purely spin-dependent gauge transformation at 
NNLO will involve only the contribution linear in G s . 
In other terms, we only need to consider the following 
transformation on the Hamiltonian 



Hamiltonian up to NNLO are then given by 
H^(r",p",S",S"*)=HW), 

i^ L °o(r",p",S",S"*) = HZoiv") {H' oN ,G s2PN }(y"), 

rW/so I II II c" ai'*\ _ zj'so f„,"\ 
-"NNLOV 1 iP 7° ) 3 ) ^NNLOV" ) 

- {-ffoN:G s 3PN} 
+ {-^olPN; G s2 Pn} 

+ {H^,G s2PN }](y"). (21) 

Note that the single prime in these equations explic- 
itly addresses the various contribution to the spin-orbit 
Hamiltonian as computed after the purely orbital canon- 
ical transformation mentioned above (note however that 
only the functional form of H^ 1PN is modified by the ac- 
tion of the orbital canonical transformation). 

Further simplifications occur in the third Poisson 
bracket of Eq. (21). First of all, since we are inter- 
ested in computing only the contribution to the spin- 
orbit interaction, the terms quadratic in spins arc ne- 
glected. In addition, from the basic relation {Si,Sj} = 
tijkSk one can show by a straightforward calculation that 
{-ff£o! G S 2pn} = (always at linear order in the spin). 
Consequently, the effect of the purely spin-dependent 
canonical transformation is fully taken into account by 
the two Poisson brackets involving the generating func- 
tions G S 2pn and G S 3pn, and the purely orbital contribu- 
tions to the Hamiltonian, H^ N and H' olPN . 

For simplicity of notation, we shall omit hereafter the 
double primes from the transformed Hamiltonian. We 
now need to connect the real Hamiltonian H to the effec- 
tive one -ffeff, which is more closely linked to the descrip- 
tion of the EOB quasigeodesic dynamics. The relation 
between the two Hamiltonians is given by [3] 
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where the real Hamiltonian H contains the rest-mass con- 
tributions Mc 2 . In terms of the nonrelativistic Hamilto- 
nian i? NR , this equation is equivalent to 
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where it is explicitly 
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H"(y") = H'(y")-{H',G s }(y"). 



(20) 



Extracting from this equation the spin-dependent terms, 
we find that the relevant terms in the new spin-orbit 



By expanding in powers of 1/c 2 up to 3PN fractional ac- 
curacy (and in powers of the spin) the exact effective 
Hamiltonian, one easily finds that the spin-orbit part 
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of the effective Hamiltonian H c g (i.e., the part which 
is linear-in-spin) reads 
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(25) 



Combining this result with the effect of the generating 
function discussed above, we get the transformed spin- 
orbit part of the effective Hamiltonian in the form as 



4ff = 4-2 (gf(S,n,p) + gf,(S*,n,p)) 



(26) 



The effective gyro-gravitomagnetic ratios g^ and g^f dif- 
fer from the ADM ones introduced above because of the 
effect of the (orbital+spin) canonical transformation and 
because of the transformation from H to H e g . They have 
the structure 



gf = 2+ ^fff NLO (a) + ±gf ™ LO (a; a, 0, 7) (27) 

gf 



if = \ + ^gf NLO (b) + ±g$?™> (b; S, C, 77), (28) 



where we made it apparent the dependence on the (in- 
dependent) NLO and NNLO gauge parameters. Includ- 
ing the new NNLO terms, they read 



gf 



a,s> 



c z 



+ aj p 2 - ( |i/ + 3a J (n • p) 2 J - -{y + a) 



9i> + — v + a + a 



(n • p) 2 ( —v v 2 + 6a - 4a - 3^ - 2 7 

1 4 16 



17 11 , 3a 

p ( -—v + y v 2 + a ~ 7 



r I" 2 + 1 + 3/3 - 37 ) p2(n • p)2 + 0§ " 2 ~ 5/3 



+ 



2 +T,p 



| + i, + 6) P 2 -(^ + 36)(„. p ) 2 -l(i + ^ + 6 



1 

+ - 

r 



1 /1 55 13 2 , 



1 59 

p~ 1 v -\ — z/" 

\4 16 2 



57 21 2 36 „, „ 
16" " ~8 ~2 + 3C " 3V ^ P (n ' P) 



7 11 i/ 2 & 

f h 7? I P" 

16 16 16 2 



15 



^ 2 -5C (n-p) 



36 
~2 



-77 + 8 - V 



(29) 



(30) 



This is the central result of the paper. The NNLO con- 
tribution to the gyro-gravitomagnetic ratios computed 
here is the crucial, new, information that it is needed to 
improve to the next PN order the spin-dependent EOB 
Hamiltonian (cither in the version of Rcf. [2] or [35] ) . Let 
us recall in this respect that in the EOB approach of [2] 
the relative dynamics can be equivalcntly represented by 
the dynamics of a spinning effective particle with effec- 
tive spin cr moving onto a ^-deformed Kerr-type metric. 
The gyro-gravitomagnetic ratios enter the definition of 



the test-spin vector cr as 

a=^(gf-2)S+ 1 -(gf.-2)S*, (31) 

that can then be inserted in Eqs. (4.16) of Rcf. [2] to get 
the spin-orbit interaction additional to the leading Kerr- 
metric part. Together with Eqs. (4.17), (4.18) and (4.19) 
of Rcf. [2] this defines the real EOB-improved, resummcd 
Hamiltonian for spinning binaries at NNLO in the spin- 
orbit interaction. 
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IV. LIMITS, CHECKS AND GAUGE FIXING 

A. The extreme-mass-ratio limit 

The effective spin-orbit Hamiltonian (26) is naturally 
connected to the test-mass (y — > 0) Hamiltonian explic- 
itly obtained 3 in [34]. To show this in a concrete case, 
let us consider the spin-orbit Hamiltonian of a spinning 
test-particle on Schwarzschild spacetime written explic- 
itly using isotropic coordinates, as given by Eq. (5.12) of 
Ref. [34]. By considering the Schwarzschild metric writ- 
ten as 



with 



ds 2 



-f(r)dt 2 + h{r)(dx 2 +dy 2 



dz 2 ) 



(32) 



where r labels here the isotropic radius 4 , r 2 = x 2 +y 2 +z 2 , 
(that is meant to be expressed in rescaled units, where 
now M ~ mi is the background mass and fi ~ mi is the 
test-particle mass), with 



„Schw 



1 



I 2r + 
1 

2r z r 



(n-p) 



16 J 



(38) 



In the v — > (Schwarzschild) limit, one has lim„_s.o {H — 
const. )/ n = lim^^o -Hcff (when dropping inessential con- 
stants), S* = Sq and S = 0. One then finds that the 
result (38) agrees in the v — > limit with Eq. (30) when 
the gauge parameters (b, S, £, 77) arc simply zero. 

In addition, in the v — > limit where the background 
is a Kerr black hole, i.e. S / 0, Eq. (29) consistently 
exhibits that both the NLO and NNLO contributions 
become pure gauge, that can just be set to zero by de- 
manding (a, a, j3, 7) to vanish. 



h= 1 



1 



2c 2 ? 



(33) 



and using rescaled variables (and making explicit the 
speed of light) Eq. (5.12) of Ref. [34] can be written as 



ISO 



— — n ISO 



(34) 



In this equation, Sq is the (rescaled) spin of the 
test-mass and we have introduced the test-mass gyro- 
gravitomagnctic ratio in isotropic coordinates , that 
is known in closed form [34] and reads 



„iso _ 
9o ~ 



where 



1 - 



2c 2 ? 



2- 



2c 2 ? 



c 2 h 



(35) 



(36) 



By transforming the Hamiltonian (34) from isotropic 
to Schwarzschild coordinates using the v — > limit of 
the (purely orbital) canonical transformation given by 
Eqs. (15)-(16), expanding in powers of 1/c 2 , (and drop- 
ping again the primes for simplicity) one obtains 



Hi 



Schw 



Schw 
,2„2 90 



n,p, S* 



(37) 



3 Note in passing that the simple procedure described in Ref. [28] 
to obtain the spin-orbit Hamiltonian is totally general and can 
be applied, in particular, to the test-mass case. 

4 Note that we use the same notation for the isotropic radius on 
Schwarzschild spacetime and the ADM radial coordinates. There 
is no ambiguity here since for the Schwarzschild spacetime ADM 
coordinates do actually coincide with isotropic coordinates 



B. Circular equatorial orbits 

Let us consider now the situation where both individ- 
ual spins are parallel (or antiparallcl) to the (rescaled) or- 
bital angular momentum vector £ = rn x p. [Note that 
in this Section the quantity r denotes the EOB radial 
coordinate (further modified by spin-dependent gauge 
terms, see below)]. In this case, circular orbits exists 
(but in the general case, when the spin vectors are not 
aligned with t, there are no circular orbits). One can 
then consistently set everywhere the radial momentum 
to zero, p r = n • p = and express the total (orbital 
plus spin-orbit part) real, PN-expanded and canonically 
transformed Hamiltonian, H(y") = H^(y") + H" so (y") 
(dropping hereafter the primes for simplicity) as a func- 
tion of r, t (using the link p 2 = £ 2 /r 2 , where £ = \£\) and 
of the two scalars a and a* measuring the projection of 
the basic spin combinations S and S* along the direction 
of the orbital angular momentum £. Following the same 
notation of [2] , we introduce here the dimensionless spin 
variables corresponding to S and S* 



cS 



cS* 
GM 2 



GM 2 ' 

and we define the projections as 

a-£ = a£, a* • t = d*£, 



(39) 



(40) 



with the scalars a and d* positive or negative depending 
on whether say a is parallel or antiparallcl to £. The se- 
quence of circular (equatorial) orbits 5 is then determined 



To avoid confusion, let us stress that we are here considering the 
circular orbits of the PN-expanded real Hamiltonian and not the 
circular orbits of the EOB-resummed real Hamiltonian, as done 
in Sec. V of Ref. [2]. This analysis is postponed to future work. 
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by the constraint 

dH(r,£,a,a*)/dr = 0, (41) 

(or equivalently by dH g/dr = 0). To start with, let 
us consider first the link between the nonrelativistic en- 
ergy (per unit mass (i) and the orbital angular momen- 
tum along circular orbits. The relevance of this quantity 
in the nonspinning case, say E c - U - C (£) = H„ (£)/fx, was 
pointed out in Ref. [38], since it provides a completely 
gauge-invariant characterization of the dynamics of cir- 
cular orbits. When the black holes are spinning, the 
same property of gauge-invariance is maintained when 
the spins are parallel (or antiparallcl) to the orbital an- 



gular momentum, so that it is meaningful to explicitly 
compute E ciic (£, a, a*) = H^ R (£)/fi + H™(£, a, a*)/n in 
this case. Since it is a gauge-invariant quantity, the result 
is independent of the canonical transformations that we 
have performed on the two-body Hamiltonian in ADM 
coordinates, so that it gives a reliable check of the pro- 
cedure we followed. As a first operation, we need to 
solve, iteratively, the constraint (41) so to obtain the 
EOB coordinate radius r in function of (£, a, a*). This 
function (that is not invariant and depends explicitly on 
the gauge parameters) reads (putting back the explicit 
double primes on r as a reminder that this is the EOB 
radial coordinate) 



r"{l,a,a*) = £ 2 { 1 + 



3 11 9 



(-9 + 3i/)l + la (33 - a(i/)) + a* - \v + b{v) 



rA 257 41 2 \ 1 1 1 / /1197 1973 3 2 „ , , 
-54 + —»- - A - + --UI — + + 6«(„) + aM + 7 (*0 



/2777 1633 7 , „. . , 

-v + —v 1 + 6b(v) + 8(v) + r}{v) 



V 16 



16 



16 



(42) 



The function E c i TC (£, a, a*) is obtained by inserting this 
relation in the expression of H{r, I, a, a*), and it reads 



1, a, a 



2£ 2 



1 1 1 



-(81-7, + , 2 ) F ---((36 + 7 ,)a+-a 



99 



2 . . 11 



cl" \\ 4 / 
-324 + 54t/- -v 2 



1701 195 



8 



(43) 



where we defined 

3861 8833 



2oi(*/) 



G4 



41 2 5 2 5 3 ,AA\ 

-k 1 v v z A v 6 , (44) 

192 32 32 64 ' v ; 



for the 3PN-accurate orbital part, with a slight abuse of 
the notation of Ref. [38]. Note that, as it should, Eq. (43) 
is totally independent of the eight gauge parameters. We 
have further verified that that the same result (43) is ob- 
tained starting from the PN-cxpandcd Hamiltonian writ- 



ten in ADM coordinates and in the center of mass frame, 
Eqs. (13)-(14). 

As a last remark, let us note that, as it was the case 
at NLO [2] , the effective gyro-gravitomagnetic ratios for 
circular orbits are gauge independent also at NNLO. To 
see this explicitly, one just imposes in Eqs. (29)- (30) the 
condition (n • p) = and the (approximate) link 



1 3 

c 2 r 2 



0(a,a* 



(45) 
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that is obtained by inverting Eq. (42) at 1PN accuracy 
and neglecting the linear-in-spin terms (that would give 
quadratic-in-spin contributions). At NNLO, one obtains 



9sl, 



9 151 



1 /51 



1 , 

r 



i 



27 
16 



39 



16 



17' (46) 



(47) 



These equations indicate that the inclusion of NNLO 
spin-orbit coupling has the effect of reducing the magni- 
tude of the gyro-gravitomagnetic ratios. The NNLO and 
NLO spin-orbit contributions act then in the same direc- 
tion, so to reduce the repulsive effect of the LO spin-orbit 
coupling which is, by itself, responsible for allowing the 
binary system to orbit on very close, and very bound, or- 
bits (see also Ref. [6] and the discussion in Sec. VI of [2]). 
We postpone to future work a detailed quantitative anal- 
ysis of the properties of the binding energy entailed by 
Eqs. (46)-(47). 



C. Gauge fixing 

We can finally exploit the flexibility introduced by 
the spin-dependent gauge transformation so to consid- 
erably simplify the expression of the effective gyro- 
gravitomagnetic ratios, Eqs. (29)-(30). This is helpful in 
the study of the dynamics of a binary system with gener- 
ically oriented spins. Reference [2] found it convenient to 



fix the NLO gauge parameters (a(v), b(v)) to 
a(u) 



3 

--v. 



u \ 5 v 



(48) 



so to suppress the dependence on p 2 at NLO. One can 
follow the same route at NNLO, i.e., by choosing the six 
gauge parameters so to suppress the terms proportional 
to p 2 , p 4 and p 2 (n-p) 2 . In this way the spin-orbit Hamil- 
tonian is expressed in a way that the circular (gauge- 
invariant) part is immediately recognizable. With (a, 6) 
fixed as per Eq. (48), one easily sees that the aforemen- 
tioned NNLO terms are removed by the following choices 
of the NNLO gauge parameters 



a{v) 

Ply) 

7(f) 
S(u) 

CM 



S- v ( 3 - v ) > 



^(13* -2), 



16 
1 

ic 
i 

ic 



(9 + 54i/-23i/ 2 ), 
(-2 + 7^ + ^ 2 ) , 



¥- 7 



(49) 
(50) 
(51) 
(52) 
(53) 
(54) 



The effective gyro-gravitomagnetic ratios are then sim- 
plified to 



9s 



2+^7 



1 5 33, 

1 / v 
~8 



1/9 3 
-r(8 + r ] 



1 



27 
16 



39 



— v 

4 



16 



-ff+ff 2 ) (n-p) 2 + ^(1 + 7.) (n-p) 4 }, 



69 
16 



57 
— i 
16 



(n-p) 



35 
16 



45 
— i 
16 



(n • p) 4 



(55) 



(56) 



This result extends the information of Eqs. (3.15a) and 
(3.15b) of Ref. [2] at NNLO accuracy. The circular- 
orbit result mentioned above is immediately recovered 
"at sight" by imposing (n-p) = 0. With this result in 
hand, one can proceed similarly to Sec. IV of Ref. [2] (as 
outlined above) to introduce the spin-dependent EOB- 
resummcd real Hamiltonian including NNLO spin-orbit 
couplings. 



V. CONCLUSIONS 

Building on the recently-computed next-to-next-to- 
lcading order PN-expanded spin-orbit Hamiltonian for 
two spinning compact objects [1], we computed the ef- 
fective gyro-gravitomagnetic ratios entering the EOB 
Hamiltonian at next-to-next to-leading order in the spin- 
orbit interaction. This result is obtained by a straightfor- 
ward extension of the procedure followed in [2] to derive 
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the NLO spin-orbit EOB Hamiltonian. We discussed in 
detail the test-particle limit and the case of equatorial 
circular orbits, when the spins are parallel or antiparallcl 
to the orbital angular momentum. In this case, one finds 
that the NNLO spin-orbit terms moderate the effect of 
the spin-orbit coupling (as the NLO terms was already 
doing [2]). 

Finally, while this paper was under review process, 
Ref. [39] appeared in the archives as a preprint: that 
study uses the Lie method to obtain effective gyro- 
gravitomagnetic coefficients that are physically equiva- 



lent to the ones presented here. In addition, it also works 
out two classes of EOB Hamiltonians that are different 
from the one considered here. 
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